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ABSTRACT 



Aims. Radiative transfer modeling of spectral lines including partial redistribution (PRD) effects requires the evaluation of the ratio 
of the emission to the absorption profile. This quantity requires a large amount of computational work if one employs the angle- 
dependent redistribution function, which prohibits its use in 3D radiative transfer computations with model atmospheres containing 
velocity fields. We aim to provide a method to compute the emission to absorption profile ratio that requires less computational work 
but retains the effect of angle-dependent scattering in the resulting line profiles. 

Methods. We present a method to compute the profile ratio that employs the angle- averaged redistribution function and wavelength 
transforms to and from the rest frame of the scattering particles. We compare the emergent line profiles of the Mg II k and Ly a lines 
computed with angle-dependent PRD, angle-averaged PRD and our new method in two representative test atmospheres. 
Results. The new method yields a good approximation of true angle-dependent profile ratio and the resulting emergent line profiles 
while keeping the computational speed and simplicity of angle- averaged PRD theory. 
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1. Introduction 

The effects of partial frequency redistribution (PRD) are im- 
portant for spectral lines that form in environments with a low 
mass density, where the average time between collisions be- 
tween atoms and electrons is long compared to the lifetime of an 
excited atomic state. Examples are the Ca II H & K, Mg II h & k 
and Lyman lines of hydrogen in the solar spectrum. Forward 
modeling of such lines requires a radiative transfer algorithm 
that includes PRD effects. 

In this paper we describe an improvement to the 
Multilevel Accelera t ed Lambda Iteration (MALI) scheme of 
IRvbicki & Hummerl (Il99llll992l) that h as been extended to in- 
clude PRD effects by lUitenbroekl (I2001I) . This improvement al- 
lows the scheme to treat PRD effects in atmospheres with strong 
velocity fields using the angle-averaged redistribution function 
without having to resort to the more accurate, but computation- 
ally expensive, angle-dependent redistribution function. 

With this new method it is possible to accurately compute 
emergent line profiles including PRD effects much faster than 
with the full angle-dependent PRD treatment. This will allow de- 
tailed comparison of radiation-MHD models with observations 
of chromospheric lines. 



2. Computation of partial redistribution effects 

The solution of the radiative transfer problem with the MALI 
scheme requires the evaluation of the ratio of the emission (ft) 
and absorption ((f)) profiles of a spectral line for which PRD ef- 



fects are important. Ignoring cross-redi stribution, this rat io in the 
observer's frame is given by (see, e. s. . IUitenbroekl200 lb : 
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Here, v and V are frequencies of the emitted and absorbed pho- 
ton in the observer's frame, n and ri their respective directions, 
y is the coherency fraction, and tij are the populations of the 
lower and upper level of the line, is the Einstein coefficient 
for radiative excitation and Pj is the total rate coefficient out of 
the upper level, Z?n is the angle-dependent observer's frame re- 
distribution function for a broadened upper and a sharp lower 
level and / is the intensity. 

The redistribution function is given by iHummerl (1 1 9621) as 
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Here g(n, ri) is the dipole scattering phase function, 
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vq)/Avd are the reduced frequencies of the ab- 



sorbed and emitted photon, with y and Avq the line center fre- 
quency and the Doppler width. The acute angle between n and ri 



1 



J. Leenaarts et al.: Approximation of angle-dependent PRD in moving atmospheres 



is given by y, H is the Voigt function and and a the line damping 
parameter. 

Note that the second term in the integral of Eq. [TJis 



t f>(v,n)I(v,n) — dv = R ij /B i 



(4) 



with Rij the upward radiative rate coefficient that appears in the 
rate equations. 

In a static atmosphere the absorption profile does not depend 
on direction. If one furthermore assumes that the radiation field 
is isotropic, i.e., 



i(y,n) = 

4-7T 



I(v, ri) dQ = J(y), 



(5) 



then the angle inte gral in Eq. Q]can be performed analytically 
(see lHubenvlll982h yielding a significantly simpler expression 
for the profile ratio: 



p(v, ri) = 1 + y 



rijBjj r 



(gn(v, v')-0(v'))/<y)dv'. 



(6) 



This simplification is often called angle-averaged PRD. 

The functio n ffn can be comput ed efficiently using the ap- 
proximation by iGouttebrozd (1 1986 1 ) with the generalization to 
cro s s -redistribution by lUitenbroekl (ll 989b . 

The assumption of an isotropic radiation field generally 
yields results very close to the full angle-dependent case 
for atmospheres w ithout, or with only a weak, velocity field 
(IUitenbroekl2002l) . 

However, Eq. [6] typically becomes inaccurate when the ve- 
locities in the atmosphere are larger than the Doppler velocity 
y/2kT/m because Eq.|5]is no longer valid, and one should ideally 
use the angle-dependent formula (Eq. [TJ. Unfortunately, using 
the angle-dependent expression is computationally expensive in 
both memory requirements and computing time. The computa- 
tion of Eq.[T]is at least a factor n a slower than the computation 
of Eq. [6] with « a the number of rays used for the direction inte- 
gration of the radiation field. 

The way the angle-dependent MALI scheme is set up one 
needs to store the intensity I(v, ri) and the profile ratio p(y, ri) in 
memory, amounting to 2 n v « a floating point numbers per spa- 
tial grid point, in addition to all other quantities that need to be 
retained. Ideally one also stores the redistribution function Ru 
to avoid recomputation every iteration, which requires an addi- 
tional rr y n\ floating point numbers per spatial grid point. 

The memory and computing time requirements are eas- 
ily met for a ID model atmosphere. However, if one wants 
to compute synthetic line profiles from time series of 3D 
radiation-MHD models, a faster algorithm to compute the angle- 
dependent PRD profile ratio is desirable, if not required. 
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Fig. 1. Plots of the radiation field used in the test computation 
as function of the reduce frequency q. Grey curves: intensity for 
different directions in the observer's frame. Black curve: angle- 
averaged radiation field in the observer's frame. Red curve: 
angle-averaged radiation field in the rest frame of a gas par- 
cel moving with a velocity (u x , u y , u z ) = (0, 1, 3) in units of the 
Doppler speed. 



the profile ratio in the rest frame, replacing the angle-averaged 
intensity in the observer's frame J(v) with the angle-averaged 
intensity in the parcel's rest frame: 

p r (v r , ri) = 1 + y^r f (sn(v r , vf) - <f>(v' t )) / f (vf) dV t . (7) 
n J p J J 



(8) 



/ r (v r ) = — (p /(v r (l - n ■ u/c), ri) dCX 



The quantity / r (v r ) can be computed incrementally from the 
observer's frame intensity I{v, ri) during the standard angle- 
frequency integration needed to compute the radiative rate co- 
efficients, without having to keep the intensity for each spatial 
location, frequency and angle in memory. 

Plugging J r (v r ) into Eq. [6] instead of J(v) yields a direction- 
independent emission to absorption ratio in the rest frame p r (v r ). 
Both J T {v r ) and p r (v r ) must be stored in memory, but as they 
are angle-independent they require a factor « a less storage than 
the angle-dependent quantities I{v, ri) and p(v, ri). The MALI 
method requires p in the observer's frame, which is accom- 
plished through a frequency shift: 

p(v, ri) = p r (v r (l + n ■ u/c)). (9) 
Numerically, this shift can be performed by a fast interpolation. 



3. Hybrid PRD 

The algorithm we propose should not be much slower than the 
angle-averaged case, yet capture the effects of velocities. It turns 
out that such a method is fairly simple: 

The agreement between angle-dependent and angle- 
averaged PRD in static atmospheres implies that the angle- 
dependence caused by anisotropy of the radiation field is minor. 
Instead, the majority of the radiation anisotropy at a given wave- 
length experienced by a parcel of gas in the observer's frame in 
a moving atmosphere is caused by sampling the line profile at 
different Doppler shifts when it receives radiation from differ- 
ent directions. If we first transform the radiation field to the rest 
frame of the moving gas parcel, we can again assume that the 
radiation field is isotropic. Then we can use Eq. [6] to compute 



4. Results 

Idealized test case To demonstrate the effect of the various 
methods to compute the profile ratio we performed an idealized 
test case computation. 

As intensity we took the emergent intensity of the Ca II H 
line computed from the standard FALC model atmosphere 
(lFontenlaetal.ll 19901) with an artificially enhanced K2R peak to 
introduce an asymmetry in the line profile. We also introduced a 
directional anisotropy by setting 
cos 6 

l(ri) = 4(1 - — ) (10) 

with I s the emergent intensity and 6 the angle between the ver- 
tical and the direction n. These intensities are shown as grey 
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curves in Fig. Q] as function of the reduced frequency q = 
(v - Vo)/Avd for different values of 9, where we assumed Avd = 
(vo/c) x 8 km s _I . The black curve shows the angle-averaged in- 
tensity in the observer's frame J. The red curve shows J r , the 
angle-averaged intensity in the rest frame of a gas parcel moving 
with a velocity of (u x , u y , u z ) = (0, 1, 3) in units of the Doppler 
spee d. The angle int egrations were performed using the A8 set 
from ICarlsonl (1 19631) . which uses 10 rays per octant. The pres- 
ence of the velocity smoothes the K2V and K2R peaks in the gas 
parcel's rest frame. The K2R peak is reduced in height and has 
an increased width, the K3 minimum is smoothed out and the 
K2V peak is a barely discernible hump at q — 4. 

Figure [2] shows p(q) given the intensity and gas velocity de- 
scribed above for two different directions (given in the upper left 
corner of the panels). We assumed y = 0.9 and njBjj/(rijPj) = 
5 x 10 5 . The black curve represents the angle-dependent case 
(AD-PRD, Eq. Q]), the blue curve the angle-averaged case (AA- 
PRD, Eq.[6j>, while the red curve shows the hybrid case (H-PRD, 
Eq. [9j. The upper panel displays scattering into the direction of 
movement and shows a blue-shifted p-profile for the AD-PRD 
and H-PRD cases. The AA-PRD case does not take the gas ve- 
locity into account and is centered at q — 0. In AD-PRD and H- 
PRD the profile ratio has a fairly smooth lowest part, reflecting 
the smoothing of the radiation field in the rest frame. Compared 
to H-PRD, the AD-PRD profile ratio exhibits additional wiggles 
owing to the directional anisotropy of the radiation. 

The lower panel shows scattering with an angle close to per- 
pendicular to the velocity. As a consequence all p profiles are 
approximately centered around p — 0. Still the H-PRD case re- 
produces the true AD-PRD profile ratio much better than the 
AA-PRD formulation. 
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: (u x ^J = (-0.79,-0.58,0.22);a= 88.60' 




1 D test atmosph eres We impleme nted the hybrid method into 
the RH code by lUitenbroekl d2001l) . which is also capable of 
using the angle-dependent and angle-independent redistribution 
function. We used this code to compute the emergent inten- 
sity from five ID plane parallel model atmospheres in non- 
LTE with the three recipes for PRD for a 4-level plus con- 
tinuum Mgll and a 5-level plus continuum HI model atom. 
As model atmospheres we took columns from a snapshot of a 
3D radiation-MHD sim ulation computed with the Bifrost code 
dGudiksen et al.ll201 The deta i ls of this particular snapshot 
can be found in iLeenaarts et al.l d2012l) . but are not important 
here: the snapshot merely provides self-consistent model atmo- 
spheres with a velocity field. No microturbulence was added. 
We treated Mgll h&k, Lya and LyyS in PRD, all other lines 
were treated using complete redistribution. For the angle inte- 
gration of the radiation field we used a ten-point Gauss-Legendre 
quadrature. 

For all atmospheres and ray directions we found that the H- 
PRD method reproduces the true AD-PRD emergent line profile 
much more accurately than the AA-PRD method. We illustrate 
this in Fig. [3] which shows the emergent intensity for a near 
vertical ray (p., - 0.95) for the three different treatments of the 
redistribution function for two different model atmospheres. 

Both atmospheres show a complex velocity and temperature 
structure (panels a and b). The second row (c and d) displays 
the emergent Mg II k intensity. For this line we expect strong 
angle-dependent PRD effects as the Doppler width of the line 
profile (2.6 km s at a temperature of 10 kK) is small compared 
to the flow velocities in the atmospheres. This is indeed the case, 
the AA-PRD profile (blue) is very different from the AD-PRD 
and H-PRD profiles. In panel c the former has a higher emission 



Fig. 2. Emission to absorption profile ratio computed using 
angle-dependent redistribution (black, Eq.[T}, angle-independent 
redistribution (blue, Eq.|6) and our intermediate method (red, see 
Sec.|2]i in the observer's frame as function of reduced frequency 
q. The cosines of the angle of the ray direction with the coor- 
dinate axes (p x ,p y ,p z ) are given in the upper left corner of the 
panel. The quantity a gives the angle between the gas velocity 
and the ray direction. 



Table 1. Comparison of running time of PRD computations. 
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0.94 
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AA-PRD 


0.073 
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0.77 
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148 


AD-PRD 


0.067 


9.40 


29.00 




491 


148 


H-PRD 


0.092 


0.06 


0.86 



peak at the blue side of the line core, whereas the latter two have 
a higher peak on the red side. In panel d the AD-PRD and H- 
PRD profiles show a stronger emission peak on the red side of 
the line core. 

The bottom row (e and f) compares the Ly a profiles. The ef- 
fect of angle-dependent PRD is much smaller than for the Mg II 
k line because the Doppler width of the line profile (13 kms 1 
at a temperature of 10 kK) is larger than the typical atmospheric 
velocities. 
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Fig. 3. Comparison of PRD treatment in two plane-parallel model atmospheres (left-hand and right hand column respectively) taken 
from a radiation-MHD model. Panels a and b show the temperature (black) and vertical velocity structure (red) of the atmospheres. 
Panels c and d show the emergent intensity in the Mg II k line, with PRD treated fully angle-dependent (black), angle-averaged 
(blue) and the hybrid approach (red). Panels e and f show a similar profile comparison, but now for the Ly a line of hydrogen. The 
stars and diamonds in panels a and b indicate the t — 1 height for the wavelengths indicated in panels c, d, e and f, with stars for the 
Mg II k line and diamonds for Ly a. 



Computational speed In Table[T]we compare the running time 
of the code for the three PRD methods for both model atoms 
in a ID atmosphere with 225 spatial points (the atmosphere in 
panel a of Fig. [3}. The quantities n v and «y RD are the number 
of frequencies used, and the number of frequencies in the PRD 
lines, respectively. The quantity fps is the time to perform the for- 
mal solution for all frequencies and angles, including the com- 
putation of J r in the H-PRD case; fpp> is the time to compute 
the profile ratios (Eqs. Q] |6]and|9]>; ?mali is the time to perform 
one full MALI ite ration, including three PRD sub-iterations (see 
lUitenbroeUiOOlh . 

The time spent in the formal solution is slightly longer in the 
H-PRD case compared to the other methods. This is due to the 
interpolations required to compute Eqs.|8]and|9]numerically. The 
computing time of the AD -PRD method is mainly spent in the 
computation of the profile ratio. The computation of the profile 
ratio is several hundred times faster in the AA-PRD and H-PRD 
cases. The total time per MALI iteration for AD-PRD is a fac- 
tor 130 (30) longer for the Mgll (HI) computation than in the 



corresponding AA-PRD computation. In contrast, the time per 
iteration for the H-PRD case is only ^10% larger than for AA- 
PRD. 

5. Discussion & conclusions 

We have presented a fast approximate method to compute the 
angle-dependent emission to absorption profile ratio needed to 
compute line profiles for which PRD effects are important. Line 
profiles computed with this method approximate the true profiles 
computed with full angle-dependent PRD very well. 

Test computations show that the new hybrid method is 30 
to 130 times faster than the angle-dependent PRD computation 
and only »T0% slower than the angle-averaged PRD method. 
This makes it possible to compute accurate line profiles from 
time-series of 3D radiation-MHD models for lines where angle- 
dependent PRD effects are important. 

Acknowledgements . JL recognizes support from the Netherlands Organization 
for Scientific Research (NWO). 
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